Exploiting environmental resonances to enhance qubit quality factors 
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We discuss dephasing times for a two-level system (including bias) coupled to a damped harmonic 
oscillator. This system is realized in measurements on solid-state Josephson qubits. It can be 
mapped to a spin-boson model with a spectral function with an approximately Lorentzian resonance. 
We diagonalize the model by means of infinitesimal unitary transformations (flow equations), and 
calculate correlation functions, dephasing rates, and qubit quality factors. We find that these depend 
strongly on the environmental resonance frequency fi; in particular, quality factors can be enhanced 
significantly by tuning to lie below the qubit frequency A. 



Introduction. - A key feature in qubit design is to gain 
good control of dephasing induced by the environment. 
A much-studied model that has yielded considerable in- 
sight into the dephasing of qubits (more generally, 2-state 
systems) is the spin-boson model Q. Most studies of 
this model assume a spectral function J{ijj) that has a 
power law form. However, several qubit systems of cur- 
rent interest are coupled to an environment which fea- 
tures rather strong resonances, which would correspond 
to a spectral function J(w) with well-defined peaks at 
characteristic frequencies. A prominent example is the 
case of flux-qubits j3| which are read out using a SQUID, 
which has a characteristic resonance frequency (of or- 
der 3 GHz), which in order of magnitude is comparable 
to the characteristic qubit energy scales (10 GHz) [^^. 

The presence of environmental resonances raises sev- 
eral interesting questions with both fundamental and 
practical implications: How is the qubit dynamics in- 
fluenced by the presence of environmental resonances? 
Can the latter be used to indirectly tune qubit proper- 
ties such as the tunneling rate or q-factor? Is it more 
advantageous to have the resonance frequency higher or 
lower than the characteristic qubit energies? 

In this Letter, we explore these questions in the frame- 
work of a model that has been used with great success to 
describe and optimize recent generations of flux qubits 
[^: it involves a spin degree of freedom (qubit) coupled 
to a harmonic oscillator with frequency O (modeling the 
environmental resonance), which in turn is coupled to a 
bath of harmonic oscillators (to provide damping) ^ . It 
can be mapped ^ onto a regular spin-boson model with 
a spectral function J(u;) featuring an almost Lorentzian 
resonance peak near Q.. We are interested not only in the 
regime where the qubit tunneling rate A is much smaller 
than Vl (which would correspond to the standard spin- 
boson model, with O playing the role of the bath cut-off 
frequency) but also in the hitherto unexplored regime 
A > il. Here standard weak-coupling, poor man scaling 
approach that predicts a downward renormalization of 
A is insufhcient; instead, we need a method sufRciently 
powerful to deal with all ratios of A/il. To this end, 
we use the flow-equation renormalization group (FER) 



approach of Wegner ''7] and of Glazek and Wilson 0. 
Interestingly, we find that A is renormalized upwards if 
the initial A is greater than 17, and that correspondingly, 
the dephasing times and q-factors are strongly increased. 
These results have the very important implication that 
hy appropriately tuning the environmental resonance fre- 
quency il, significant additional control of the qubit dy- 
namics can indeed by obtained. 

Spin-boson-model. - We consider the Hamiltonian 



H = -ja^ + ^a, + {B^ +B) 



kV"k 



nB^B + Y^^J>lb, + iB^ + BfY.^^ 



(1) 



which describes a 2-state-system with asymmetry energy 
e and tunneling matrix element A, coupled linearly with 
strength 5 to a harmonic oscillator with frequency fl, 
which is itself linearly coupled with strengths Kk to a bath 
of harmonic oscillators. The coupling to the environment 
is completely defined by the spectral function J(lu) = 
X^k ^k'^(^ ~ i^k) = Tu!Q{u!c — u>), which is as usual taken 
to be of ohmic form to model the dissipativc environment. 
This s_Y_stem can be mapped to a spin-boson model of the 
form 



n 



1 'sr^ 

— cr, 
2 

k 



Ak(4 + \)+E^k4^k. (2) 



where spin dynamics depends only on the structured 
spectral function J{uj) = J^k ^k'^(^ ^ ^k) given by 



with a = 



8iy 



(3) 



Flow equation renormalization.- Flow equation renor- 
malization (FER) is based on infinitesimal unitary trans- 
formations of the Hamiltonian We follow the ap- 
proach of Ref. 0, and mention only the main steps here: 
(a) In order to decouple the two-level system from 
its environment we apply a sequence of unitary trans- 
formations U{1) to Eq.©: n{l) = U{l)nW{l). Here 
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7i(Z = 0) = 7i is the initial Hamiltonian, 'H{1 = oo) is 
the final, diagonal Hamiltonian and I denotes the flow 
parameter, which characterizes the square of the inverse 
energy scale being decoupled. In differential formulation 
this transformation reads 



dnjl) 
dl 



= [77(0,H(0] with 



Vil)-'^U-\1). (4) 



(b) The canonical choice for the generator rj suggested 
by Wegner is rjc ~ [WojW] with Hq = --jcTx + §CTz + 
Sk^k^k^k However, since r\c generates new coupling 
terms (originally not present in the Hamiltonian) it is ad- 
visable to modify our generator. For e 7^ we choose ^3 



k k 

+ E '?^(^ - ^I) + E ^kq(&k + bl)ib, - bi), 

k kq 

for which Eq.l^J closes for terms linear in bosonic 
operators. We neglect small higher-order terms in 
[ry, Ti.] that contain a coupling of the system to two 
bosonic modes. The parameters ?7k and T^kq in Eq.© 
are given by rj^ ^ -(Ak/2)(eA/wk)/(t^k, 0> '7k = 
-(Ak/2)A/K,0, = ~iXU2)[iiol-e')/Lo^]fiu;^J), 

and ?7kq - AV(2Ae)tanh(/3Ae/2)Ak_AqWq/(wk " 
Lul)[f{uj]^, I) + f{LUq, I)] with Ae \/A2 +e^. We choose 

/(wk, = [w^K - A,)] / [A2(u;k + A,)] . By compar- 
ing numerical results for the £ = (see 10]) and the £ 7^ 
Ansatz we see that for £ 7^ 0, due to our particular choice 
of /(wk, I), we are restricted to couplings a < 0.02 which 
is a reasonable bound for experimental realizations. For 
an alternative Ansatz (for £ 7^ 0) see also Ref. \vi\ . 

(c) Eqs.l^l and (0) give us a set of differential equations 
(flow equations) for the parameters in the Hamiltonian, 
namely e{l), A{1), and Ak(0 [respectively J{ujJ)]: 

-9/ A/A = Jduj coth ^ J{uj,l)f{uj,l), 9z£ = 0, (6) 

diJ{uj,l) = -2 f {uj, I) {lo^ ^ A/) J {lj, I) + tanh^ X 
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J{uj,l) / du! 



,uj'J{uj',l) 



[f{u;,l) + f{cu',l)]. (7) 



Note that according to 10) the bias e is not renormalized. 
The renormalization of the bath frequencies Wk vanishes 
in the thermodynamic limit of infinitely many modes. 

(d) Observables such as have to be subject to the 
same sequence of infinitesimal transformations as the 
Hamiltonian: daz{l)/dl = [ri{l),az{l)]- For the flow of 
CTz we make the Ansatz 

a^l) = hil)a. + sil)a^ + r{l) + icTyJ2filil)ib^-bl) 

k 

+ J2^^.xui)+^^xmiK+bi) (8) 
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FIG. 1: Dephasing times for an ohmic bath with spectral 
function J(aj) = 2aujQ{u!c — a;) as a function of a. The PER 
result [e = and uic = lOA(O)] is compared with results from 
RTRG calculations El and WCC _14:]. The inset shows a 
typical FER spin-spin correlation function. 



We neglect (small) terms in [77, a^] which contain a cou- 
pling to two bosonic modes. The calculation of the flow 
equations for the six parameters h, s, r, xjj, and /i^ 
in Eq.® is straightforward (31 . 

(e) To calculate correlation functions of the form 
C{t) = 5(o'z(i)crz(0) -f crz(0)(Tz(i)) we use the decou- 
pled Hamiltonian 7i(cx)) = — [A(oo)/2]crx + [e/2]crz -I- 
Ek ^k^k^k' and Eq.® for I = 00. For T = the Fourier 
transform C{u}) of C{t) then takes the form (here all pa- 
rameters are taken at I = 00): ^3 
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(5(CJ - CJk). 



(9) 



Numerically, one finds that h(oo) = s{oo) = and i) 
r{oo) ~ for e = or ii) r(oo) 7^ for e 7^ 0. Therefore, 
of the terms in the first line of © only 6{oj) remains and 
describes the nonzero expectation value of for systems 
with asymmetry. 

In order to obtain quantitative results for the corre- 
lation function C{uj), we numerically integrate the fiow 
equations up to some value Zq, which is taken sufficiently 
large that the final results do not depend on it. From the 
numerical results for C{uj), which reflects the dynamics 
of the two-level system, we extract "dephasing times", 
defined as the widths at half maximum of the resonances 
occurring in C(w) (as depicted in the inset of Fig.QJ. For 
zero bias (e = 0) a sum rule of the form dujC{uj) = 1 
should hold p. 

Ohmic bath: We start by comparing our results for the 
dephasing time for a spin-boson model with an ohmic 
bath, J{lij) = 2aujQ{L0c — lo), with results from real-time 
renormalization group (RTRG) |l3l | and weak coupling 
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FIG. 2: Spin-spin correlation function as a function of fre- 
quency for experimentally relevant parameters discussed in 
Ref. II: a = 0.0006, A(0) = 4 GHz, e = (this is the so- 
called "idle state"), f7 = 3 GHz, V = 0.02, and = 8 GHz. 
The sum rule is fulfilled with an error of less than 1%. (a) 
Blow up of the peak region reveals a double peak; (b) blow 
up of the larger peak, (c) term scheme of a two level sys- 
tem coupled to an harmonic oscillator, drawn for A(0) ^ f2; 
(a = 0.0006 corresponds to g/fi ^ 0.06.) 

calculations (WCC) Q- Figure n shows the dephasing 
time T as a function of the coupling strength a. For weak 
coupling the dephasing time (at T = 0) is given by |Q| 



4/J[Ae(oo)]. 



(10) 



We find very good agreement with RTRG and WCC. 

Structured bath/weak coupling: We now turn to the 
structured spectral density given by Eq.((31). The main 
features of the corresponding system [Eq.Q])] can al- 
ready be understood by analyzing only the coupled two- 
level-harmonic oscillator system (without damping, i.e. 
r = 0). For e = this system exhibits two charac- 
teristic frequencies, close to ft and A, associated with 
the transitions 1 and 2 in Fig. [Sfc). These should also 
show up in the correlation function C(w); and indeed 
Fig. Efa) displays a double-peak structure with the peak 
separation somewhat larger than (A — fi), due to level 
repulsion. The coupling to the bath will in general lead 
to a broadening of the resonances and an enhancement 
of the repulsion of the two energies. Due to the very 
small coupling (a — 0.0006) peak positions of C{uj) in 
Fig.Elcan with very good accuracy be derived from a sec- 
ond order perturbation calculation for the coupled two- 
level-harmonic oscillator system, yielding the following 
transition frequencies [depicted in inset (c) of Fig. [2]: 
uji,+ - wo,+ = n- g22A(0)/[A2(0) - n^] w 0.98717 and 
Wo - - uja.+ = A(0) + .g22A(0)/[A2(0) - n^] w 1.3461). 
With the two peaks we associate two different dephas- 
ing times, TQ and ta, as shown in inset (a) and (b) of 
Fig. 13 In Fig.|2Ia) both these dephasing times are shown 
as functions of A(0)/r2 for a = 0.0006. Moreover, ta is 
compared to the WCC result Tw of Ea. (fTn|l . This com- 
parison is expected to work well for J(17)/0 <C 1, and 
indeed it does [here J{n)/n « 0.08]. Due to the small 
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FIG. 3: FER results for dephasing times (ta and tq) for the 
structured bath [Eq.||3J] compared to results from WCC (tw) 
given by Eg. llOL with A(cxd) occurring therein calculated us- 
ing FER. Parameters are e = 0, = 3 GHz, F = 0.02. (a) 
weak coupling: a = 0.0006 and lJc = 8 GHz. Sum rules are 
fulfilled with errors of less than 1%. Inset: Spin-spin cor- 
relation function for A(0) = 0.98712. (b) stronger coupling: 
a — 0.01 and Uc = 9 GHz. Sum rules are fulfilled with errors 
of less than 3%. Inset: x - Renormalized tunneling matrix 
element; o - quality factor q — rAA(cxD)/2. 

coupling between two-level-system and harmonic oscill- 
tor, g « 0.0611, the dependence of ro on A(0) is very 
weak. 

For A(0) close to 17 the two resonances in C{u!) merge 
to a symmetric double-peak structure as shown in the in- 
set of Fig.l^^a). Here a characterization by two different 
time scales becomes difficult. Therefore the correspond- 
ing data points in Fig. Ola) have not been included. 

Stronger coupling to bath: Figure |3Ib) shows ta, tq, 
and Tw for a larger coupling strength of a = 0.01. Fig- 
ure EI a) shows one of the calculated correlation func- 
tions. Note that the stronger coupling a leads to a larger 
separation, or "level repulsion", between the A- and 17- 
peaks than in Fig. |21 The inset of Fig. Olb) shows the 
renormalized tunneling matrix element A(c») as a func- 
tion the initial matrix element A(0). Very importantly, 
for A(0) ^ 17, A increases during the flow, whereas for 
A(0) ^ 17, it decreases This behavior can be under- 
stood from the fact that /(w, I) in Eq.® changes sign at 
uj = A: If the weight of J{w) under the integral in © 
is larger for uj > A(0), which is the case if A < 17, then 
positive values of f{w,l) dominate and A{1) decreases 
[A(cx)) < A(0)]; conversely if the weight is larger for 
u; < A(0), A{1) increases [A(oo) > A(0)]. Note also, 
that the upward renormalization towards larger A(oo) 
in the inset of Fig. ISJb) is stronger than the downward 
one towards smaller values, i.e., the renormalization is 
not symmetric with respect to A(0) = 17. The reason 
for this asymmetry lies in the fact that /(w, I) has a 
larger weight for w < A than for cu > A. Also ta and 
even Tw ~ 1/ J[A(oo)] in Fig. IJl^b) show an asymmetric 
behavior with a steep increase at A(0) « 17: dephas- 
ing times for A(0) > 17 are larger than for A(0) < 17. 
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FIG. 4: Spin-spin correlation function for the structured bath 
[Eq.JSJ] as a function of frequency [iq . The maximum height 
of the middle peak in (b) is ~ 7.2. 

That this happens, although J is more or less symmet- 
ric around its maximum, is a direct consequence of the 
stronger renormalization of A for the latter case. Also 
the quality factor [g-factor, see inset in Fig.lS^b)] defined 
as q — taA(cx))/2 shows this asymmetric behavior (be- 
ing larger for A(0) > 57 than for A(0) < O) with a steep 
increase at A(0) = ft from 8 to 43, i.e., by a remark- 
ably large factor of w 5. We consider the asymmetry of 
the renormalized tunneling matrix element, the dephas- 
ing time and the quality factor as the central results of 
this Letter: By tuning H. such that A(0) > fl, dephas- 
ing times can be significantly enhanced (as compared to 
A(0) < n). Note also that to in Fig. I^l^b) shows a much 
stronger dependence on A(0) than in Fig. Ol^a). This is 
due to the stronger coupling {g w 0.3f2) of the two-level 
system to the harmonic oscillator in (b). 

Nonzero bias: We now turn to the case of nonzero bias, 
e 7^ 0. A second order perturbation calculation, analo- 
gous to the zero bias case, shows that a third resonance 
in C{u!) is expected to show up at an energy scale A^ + Q. 
Indeed it does, as exemplified in FigureHfb ), w hich shows 
a typical result for C{uj) for nonzero bias jl^- With ev- 
ery resonance we associate a dephasing time (analogous 
to the zero bias case). Figure l^a) shows all three de- 
phasing times (ta, tq, and ta+q) as a function of e. ta 
is compared to the weak coupling result r^. As expected, 
Ta shows a minimum at ~ -^/fi^ — A^{Q), which cor- 
responds to the maximum of J{uj,l — 0). Beyond Emin, 
TA increases whereas tq and r^+n decrease. This is the 
expected behavior: In the limit e ^ oo (A ^ 0), C(t) 
respectively C{uj) should become independent of all bath 
characteristics, i.e., C{t) 1 and C{llj) — > S{u!). In this 
limit the dephasing times should show the following be- 
havior: TA oo. To ^ and TA+n 0. In Fig. |SJb) 
and (c) the renormalized tunneling matrix element and 
the quality factor are shown as function of the bias. Note 
that since e is not renormalized [see Eq.®], Ae(oo) as 
a function of e/fl docs not show a strong asymmetry, in 
contrast to the case e = 0. As a direct consequence, de- 
phasing times and quality factors do not change much at 
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FIG. 5: (a) FER results for dephasing times (ta, tci, and 
TA+n) for the structured bath [Eq.Q], compared to results 
from WCC (tw) given by Eg. llOt . (b) Renormalized tunneling 
matrix element, (c) q-factor q — TAAe(oo)/2. Parameters: 
a = 0.01, A(0) = 2 GHz, Q, = 2.75 GHz, F = 0.02, and 
ujc^9 GHz. 



Eniin- Finally, t and g as a function of A(0) for fixed e can 
be shown |l2| to show a qualitatively similar behavior to 
Fig. El 

Summary - We used flow equation renormalization to 
study a 2-level-system coupled to a damped harmonic 
oscillator for arbitrary ratios of A/il. We find that by 
tuning the system into the regime A > il, which is stud- 
ied here for the first time, dephasing times and q-factors 
can be significantly enhanced. 
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